Group delay of electromagnetic pulses through multilayer 
dielectric mirrors beyond special relativity 
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■ A number of unexpected phenomena historically broke through corre- 
sponding the traditional physics framework, thus motivating scientists to 
develop more advanced and accurate theories. A well-known example is 

%— > \ 

Q-i quantum tunnelling in a-decay. The quantum tunnelling rate can be de- 

O. 

c/3 ■ scribed accurately by the quantum theory but cannot be explained by clas- 

o: 

sical mechanics. In quantum tunnelling, the length of time during which a 
particle tunnels through a barrier remains undetermined®^. Based on clas- 
i— —>. sical quantum and classical electrodynamics theories, the group delay during 

quantum tunnelling is independent of barrier thickness. However, in more 

>: 

accurate theories, such as the general theory of relativity 7 , is this condition 
still valid? In this Letter, we investigate the group delay of optical pulses 

■ through multilayer dielectric mirrors (MDM) combined with gravitational 
wave (GW)I™1 We find that the delay increases linearly with MDM length 

l ■ 

for the transmitted wave packet, whereas the Hartman effect disappears. 

^ ■ 

Our study provides insight into the nature of both quantum tunnelling and 
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GW. 

The length of time during which quantum particles tunnel through a barrier has 
attracted considerable attention for both fundamental and technological reasons since 
the 1930s^^. Hartman calculated the tunneling of a wavepacket through a rectangular 
potential barrier^ and found that group delay becomes constant as barrier length in- 
creases. This phenomenon, known as the Hartman effect, implies that for sufficiently 
large barriers, the effective group velocity of a particle may be superluminal. Although a 
number of experiments have reported observations of electromagnetic waves propagating 
with "superluminal tunneling velocities," the definition of tunneling time and its exact 
physical meaning based on experimental results remain under heated debate^*'. A large 
number of tunneling time definitions have been proposed, including group delay or phase 
time^, dwell time^, Larmor times^n^, and Biittiker-Landauer time^. Winful recently 
proposed that the group delay in tunneling represents a lifetime of stored energy escap- 
ing through both sides of the barrier and does not represent a transit tim o^^l Thus, 
the issues of superluminality, causality, or the speed of information transfer do not even 
arise. 

Beyond Schrodinger's nonrelativistic quantum mechanics, the group delay for Dirac 
particles traveling through a potential well was also studied by using Dirac's fully rela- 
tivistic quantum theory^. The behavior of Dirac particles is found to be the same as 
that in nonrelativistic quantum mechanics. Liu et al. recently studied the effect of the 
electromagnetic fields of the group delay of electrons and found that the group delay 
of the transmitted wave packet increases linearly with barrier length for the transmit- 
ted wave packelP^. This peculiar tunneling effect is attributed to current leakage in a 
time-dependent barrier generated via the electromagnetic fields^"^. If the quantum 
fluctuation or the zero-point field is considered, all potential barriers are combined with 
electromagnetic fields. Thus, in the framework of quantum field theory, the Hartman 
effect of electron tunneling may disappear because of the existence of electromagnetic 
vacuum fields. 

However, in photon tunnelling, a number of theories based on the special relativis- 
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tic covariant Maxwell equations have proven that the group delay of photons becomes 
constant as the length of optical structures incxeases^™^"^. However, whether the 
Hartman effect still exists when a more accurate theory (e.g. general relativity theory) 
remains unclear. One of the unique predictions in general relativity theory is the exis- 
tence of gravitational wave (GW)PH^. If optical structures such as a multilayer dielectric 
mirror (MDM) (i.e. a one-dimensional photonic crystal) is irradiated by GW, the center 
frequency and the width of the photonic bandgap of the MDM will vary with the GW. 
Similar to the electron tunnelling in a time-dependent barrier, variations of the photonic 
bandgap will result in an additional leak photon current. Such current may propagate 
at the speed of light. Thus, the Hartman effect may also be absent in photon tunnelling 
within the framework of general relativity theory. 

In this paper, we investigate the effect of GW on the group delay of photon passing 
through MDM. Our simulation shows that with a thick MDM, group delay increases 
linearly with increasing barrier width. The group velocity is approximately 2.95 x 10 8 
m/s, slightly less than the speed of light in in vacuum. Superluminality or causality no 
longer occurs. We also find that the group delay of tunneling photons is sensitive to GW. 
Our study may facilitate further understanding of both quantum tunnelling and GW as 
well as provide a different method for the detection of GW. In particularly, the MDM 
comprises alternating dielectric layers and vacuum layers [see Fig. 1(a)]. All layers are 
nonmagnetic ([/, = 1), and the thicknesses of the dielectric layers (vacuum layers) satisfy 
D\ = Ao/4y / ?i (L>2 = A0/4 + C2A0/2), where e\ is the permittivity of dielectric layers, Ao 
is the center frequency of the input electromagnetic pulse, and C2 is a positive integer. 
The group delay of tunneling photons is generally more sensitive to GW at large £2- 
We set C2 = 5 in this paper, unless otherwise specified. The electromagnetic pulse is 
incident along the normal of the surface of MDM, and the propagation direction of the 
plane-polarized GW is parallel to the surface of MDM. 

When the polarization GW occurs, the layer spacing of dielectric layers will vary 
with the GW. Similar to the case of two masses separated by a distance D along the Z 
direction that are coupled by a lossy spring, the equation for the differential motion of 
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the masses in the GW becomes^"2lH2Z] 

d 2 z R uj m0 dz R 2 1 ^^22 n M s 

where zr is the proper relative displacement of the two masses, h 22 is the the perturbation 
matrix (tensor) element resulting from the GW; and ui m Q and Q are the natural frequency 
and the q-factor of the lossy spring, respectively. If the dielectric layers in MDM are 
free and under zero-velocity initial conditions, the equation of differential motion can be 
described by™22 

zr/D = h 22 /2 = A GW cos(u;Gwt)/2, (2) 

where Aqw is the GW amplitude. Thus, in the propagation of electromagnetic waves, 
the permittivity distribution will also change with time. To study such a time-dependent 
photon scattering process, we employ the finite-difference time-domain (FDTD) method 
to solve the time-dependent Maxwell equations numerically^"^. FDTD is a time- 
domain technique that can provide animated displays of electromagnetic field movement 
through various models, e.g., the photonic crystals, radar, etc. In the FDTD method, 
the one-dimensional Maxwell equations are replaced by a finite set of finite differential 
equations^ 11 ^ 1 



E^(k) = E^(k) - \H n y +1/ \k + 1/2) - H n y +l/ \k - 1/2) 

H^ +1/ \k + 1/2) = H^ 1/2 (k + 1/2) - 7IK £L_ [E n + i [k + 1} _ E n {k)] 

where E x (H y ) is the electric field (magnetic field) of the electromagnetic wave, (i, k) = 
(i Ax, kAt) denote a grid point of the space and time; and for any function of space 
and time F(iAx,kAt) = F k (i), Tqw depends on the GW wave-induced displacement. 
According to Eq. [2j in the vacuum layers Tqw = 1 + Aqw cos(uiQwt) /2. The thick- 
ness changes of the dielectric layers are minimal because the natural frequency of di- 
electric layer is nonresonant with the frequency of GW, and the thickness the dielec- 
tric layers is considerably less than that of the vacuum layers with large C 2 - Thus, 
we set Tqw = 1 in dielectric layers. At the input boundary, a Gaussian electro- 
magnetic wave packet is injected. The wave function at the input boundary is set as 
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E x = H y = exp — Air (t — To) 2 /tq e tuot . To reduce distortion and numerical errors, 
a relatively long pulse is used: tq = 200To, where To is the period of the electromagnetic 
wave. 

By numerically solving Eq. Q directly in the time domain, the propagation of a 
wave packet through a barrier can be demonstrated in real time. For computational 
stability, the space increment Ax and the time increment At need to satisfy the relation 
Ax > cAt. Furthermore, the space increment Ax must be significantly smaller than 
the wavelength of electromagnetic wave Ax < Aq/12, and the time increment At must 
be considerably smaller than the period of the GW. To ensure high precision, the space 
increment Ax = Ao/1.5 x 10 3 and the time increment At = 2 x 10~ 4 To are used. When 
the space and time increments are increased or reduced 10 times, the error is less than 
3%. 

Numerical results of the group delay (i.e., the delays of the peaks of the reflected and 
transmitted pulses) are shown in Figs. l(b)-(f). Fig. 1(b) shows the group delay for the 
reflected wave packet as a function of the number of MDM periods umdm- Similar to 
the traditional quantum tunneling, the group delay is saturated by increasing umdm, 
and the saturated group delay is identical to the dwell time. Meanwhile, the group delay 
for the reflected wave packet is unaffected by the extrinsic GW. 

However, the influence of the GW on the group delay of a transmitted wave packet 
differs from the given case. As shown in Fig. 1(c), for the GW amplitude Aqw = 1 x 
10~ 4 , the group delay for the transmitted wave packet increases linearly when umdm > 
24. This result can be explained by the variations of the photonic bandgap attributed 
to the GW. Similar to the electron tunnelling in a time-dependent barrieif^^, the 
variations of the photonic bandgap will result in an additional leak photon current. Given 
that the amplitude of the additional leakage current attributed GW is quite small, for 
the case of a small number of MDM periods, e.g., umdm < 15, the tunneling current is 
significantly larger than the additional leakage current and is unaffected by the extrinsic 
GW [see Fig. 1(c)]. However, for the case of large number of MDM periods, e.g., 
n MDM > 24, the tunneling current is significantly weaker than the additional leakage 
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Figure 1: (a) Schematic diagram of the tunneling process in an MDM structure, (b) 
Group delay for the reflected wave packet and (c) that for the transmitted wave packet 
as a function of the number of MDM periods with different GW amplitudes. Incident 
(black lines), tunneled (red lines) and reflected (green lines) pulses with GW amplitude 
Agw = 1 x 10 -4 for the following number of MDM periods: (d)nMDM = 3, (ejriMDM = 
19, and (f) umdm = 30. The inset shows the normalized tunneled pulse overlaid with 
the incident pulse. 

current. The additional leakage current determines the group delay, and the group delay 
increases with the amplitude of GW. Under a weaker GW, the additional leakage current 
will determine the group delay for the transmitted wave packet with a relatively larger 
nMDM- Specifically, for Aqw = 1 x 10~ 4 (Aqw = 2 x 1CT 5 ), the group delay increases 
linearly when umdm > 24 (umdm > 26). Thus, if umdm is sufficiently large, even under 
a relatively weak GW (e.g., the GW background radiation^}), no Hartman effect occurs. 
On the other hand, for a non-strictly periodic GW emitted by various sources(e.g., the 
chaos compact binary systeirPS), the time-dependent variations attributed to the GW 
will also result in an additional leak photon current, thus modifying the group delay. 
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From Fig. 1(c), we can also find that the group velocity of the additional leakage current 
is independent of the amplitude of GW. The group velocity is approximately 2.95 x 10 8 
m/s, which is the same as the speed of light in vacuum. No superluminal appears. 

The transmitted wave packet includes two components: the tunneling current and 
the additional leakage current induced by the time-dependent modulation, which may 
distort the transmitted wave packet. However, the relative weight of the contribution of 
the tunneling current and that of the additional leakage current vary with the number 
of MDM periods umdm- For a small number of MDM periods, e.g., umdm = 3, the 
amplitude of the tunneling current is significantly larger than that of the additional 
leakage current. Similar to in traditional quantum tunneling, the distortion is minimal 
[see Fig. 1(d)]. As the number of MDM periods increases, the tunneling rates decrease 
rapidly. Meanwhile, the number of MDM periods only has a slight effect on the additional 
leakage current. For a large number of MDM periods, e.g, umdm = 19, the amplitude 
of the additional leakage current and that of the tunneling current are comparable. 
Considering that the group delay of the additional leakage current is larger than that of 
the tunneling current, a serious distortion of transmitted wave packet occurs [see Fig. 
3(e)]. However, if the number of MDM periods is sufficient, e.g., umdm = 30, the 
tunneling rate is very small, and the additional leakage current is the main contributor 
to the transmitted wave packet. No distortion occurs at this scale [see Fig. 3(f)]. 

However, we have to determine whether the group delay of the undistorted trans- 
mitted wave packet in a thick MDM with GW equates to tunneling time. In traditional 
quantum tunneling, the consensus is that the group delay does not equate to a tunneling 
time for the following reasons: the group delay is equal to the dwell time, which denotes 
a lifetime of stored energy escaping through the barrier, not a pulse traveling through 
the barrier^"^; to reduce the distortion, the length of the wave packet is often signif- 
icantly larger than the gap length of the barrier^M^ZI, and the peak of the input pulse 
does not even propagate into the barrier. However, the group delay of photons with GW 
is different. The group delay with GW cannot be explained by the dwell time because 
such delay is considerably larger than the dwell time. The full width at half maximum 
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Figure 2: (a) Group delay of the transmitted wave packet as a function of the number 
of MDM periods for different GW frequencies with Aqw = 1 x 1CP 4 . (b) Group delay of 
the transmitted wave packet as a function of the number of MDM periods for different 
vacuum layer thicknesses with Aqw = 1 x 10 -4 , ojqw = O.Iloq, and To = 2000To. The 
inset shows the relative group delay as a function of the number of MDM periods. 

(FWHM) of the injected wave packet is approximately 66 Ao (or 66To), which is smaller 
than the total optical path length of MDM at umdm = 30. Meanwhile, the group delay 
for umdm = 30 is approximately 80To, which is larger than the FWHM of the injected 
wave packet. The peaks of the injected wave packet and that of the transmitted wave 
packet are distinguishable. Thus, the group delay of the transmitted wave packet in 
a thick MDM with GW may be regarded as the tunneling time. However, although 
the group delay equates to the tunneling time, the tunneling time is increased with the 
number of MDM periods. Thus, the super luminality in the tunneling process no longer 
occurs. 

Notably, the group delay in the tunneling process also shows good sensitivity to the 
GW. For instance, for umdm = 25 and Aqw = 1 x 10~ 4 , the group delay without (with) 
GW is approximately 5.3To (65To). The group delay is increased by approximately 12 
times. In the GW detection by using a Michelson interferometer, the interferometer 
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measures the intensity rather than the time delay of the interference light. We can 
still make a comparison with Michelson interferometer. For a Michelson interferometer 
with an arm length of L Q = 75Ao (same as the total optical path length of MDM with 
nMDM = 25), only when the GW amplitude Aqw is approximately 1.6 x 10~ 3 (i.e., the 
GW amplitude satisfies 4irAcwL o /^0 ~ 7r /2), the intensity of the interference light with 
and without GW can vary by approximately 12 times. 

The sensitivity of the group delay to the GW depends on the frequency of the GW 
[see Fig. 2(a)]. For high-frequency GW, e.g., ojgw = ^o> the variations of the photonic 
bandgap occur rapidly, the stored electromagnetic energy in MDM changes slightly, 
and the additional leak photon current is relatively weak. The group delay increases 
remarkably only when hmdm > 22. However, for a relatively low-frequency GW, the 
additional leak photon current is enhanced. The group delay increases remarkably when 
nMDM > 21 (hmdm > 19) for ojqw = 0.25^o (^GW = O.lwo). For an extremely low- 
frequency GW, e.g., ujgw = O.Olwo, the period of the GW becomes larger than the 
FWHM of electromagnetic waves, and the effect of the GW on the group delay becomes 
small. For umdm < 25, the group delay does not change significantly. 

The sensitivity also depends on the thickness of the vacuum layers. For MDM with 
thicker vacuum layers (i.e., large £2), the GW-induced variation of layer spacing is en- 
hanced, and a larger additional leak photon current can be achieved. As shown in Fig. 
2(b), group delay increases remarkably when hmdm > 12 for (2 = 25. However, for 
C2 = 5, group delay increases remarkably only when hmdm > 19. On the other hand, 
the relative group delay tdt/L\ of the additional leakage current is independent of C2, 
where L\ Q is the optical path length of each MDM period [see the inset of Fig. 2(b)], 
which indicates that the pulses propagate with the same group velocity for different £2- 

Finally, we discuss the experimental realization of our theoretical predication. Al- 
though a strong GW is used in the numerical calculation, our results show that the group 
delay of the tunneling photons is sensitive to GW and that sensitivity can be increased 
for a relatively low- frequency GW or for thick vacuum layers. Owing to ultrafast laser 
technology, an extremely high time resolution can be achieved. Thus, the detection of 
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the effect of the GW on group delay may be feasible. Notably, a pulsar emits a beam 
of electromagnetic radiation with ultra-high accuracy and stability. Thus, when the 
pulsar electromagnetic radiation tunnels through a strong GW radiation source such as 
the black hole binary, the variation of the group delay should be detectable through 
astronomical observation. 

In conclusion, we have calculated the group delay of optical pulses through MDM 
combined with GW. We found that the group delay increases linearly with MDM length 
for the transmitted wave packet. The Hartman effect disappears. This peculiar tunnel- 
ing effect is attributed to the additional current leakage attributed to the GW-induced 
variations of the photonic bandgap. We also show that the group delay of the tunneling 
photons is sensitive to GW. For a relatively low-frequency GW or thick vacuum layers, 
the sensitivity can be enhanced remarkably. Our study provides insight into the nature 
both of the quantum tunnelling and GW as well as a novel process by which to detect 
the GW. 

This work was supported by the NSFC Grant Nos. 10904059, 11004199, 11104232, 
11173012, and 11264030, the NSF from the Jiangxi Province Nos. 20122BAB212003. 
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